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Abstract 

We discuss optical forces exerted on particles, either dielectric or metallic, near a subwavelength 
slit illuminated by a photonic nanojet. We compare those cases in which the Mie resonances are 
or are not excited. The configurations on study are 2D, hence those particles are infinite cylinders 
and, in order to obtain extraordinary transmission, the illuminating beam is p-polarized. We show 
the different effects of these particle resonances on the optical forces: while whispering gallery 
modes under those illumination conditions weaken the force strength, this latter is enhanced by 
localized plasmon excitation. Also, illuminating the slit with a nanojet enhances the optical forces 
on the particle at the exit of the aperture by a factor between 3 and 10 compared with illumination 
of the slit with a Gausian beam. In addition, the pulling force that such a small resonant metallic 
particle suffers on direct illumination by a nanojet, can change by the presence of the slit, so that 
it may become repulsive at certain lateral positions of the particle. 
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INTRODUCTION 



Studies on optical forces on micro and nano-objects, both in their apphcations to trapping 
[l-iol as weU as to optical binding 11-13|, show the sensitivity of these techniques to the 
thermal action on the kinetics of these systems. This involves to increase both the numerical 

n 

aperture NA and the power of the illuminating beam [9'] to control the experiments; the 
former procedure has obvious limitations for the manipulation of biological specimens. 

The effectiveness of optical trapping increases 0, Q by illuminating particles through 
subwavelength apertures in supertransmission (l5l-l26l|. namely on excitation of the aperture 
morphology dependent resonances (MDR). This allows lower illumination powers and its 



performance is greatly improved 



14l | when also the Mie resonances of the particles are 



excited, i.e. their whispering gallery modes (WGM) 

(LSP) |33l-l36| . This enhances the aperture transmittance and localization of light j37 
Nevertheless, illuminating with photonic nanojets (PNJ) j40l-l43| constitutes an alterna- 



27H32| or localized surface plasmons 



tive means to enhance that transmission and 
PNJs have subwavelength spatial resolution 



oca 



41 



ization of light through apertures |44l |. 



461] and hence are of great interest 



for microscopy and detection at the nanoscale [47"]. Since however they are nonresonant fo- 
cusing effects, their appearence is not so narrowly constrained by the constitutive parameters 
and morphology of the particles as Mie resonances are. 

As far as we know, in contrast with the extensive study on nanojets, there is only one 
report concerning optical forces on a metallic particle close to a PNJ j48|. We thus present 
here an analysis on the effect that the presence of subwavelength apertures has on these 
forces; a subject that combines the phenomenon of extraordinary transmission with the 
nanojet focusing and localization of optical energy, and which has never been addressed. 
The nanoparticle is either dielectric or metallic, and its MDRs may be excited. 

We shall employ 2D configurations, so that the particles are cylinders with cross section 
in the plane of calculations, and the apertures are slits. It is well-known that this accounts 
for the main features of the phenomenon providing one does not look for depolarization 
effects l49l-l5l|and it is certainly the case for PNJs j40l-l43| . Therefore in order to obtain 



slit supertransmission, the incident wave is p-polarized. 



We will employ the Maxwell stress tensor (MST) 



BB,[i4, 



50I-I55I] . The fields are calculated 



by a finite element method p^]. Section 2 accounts for these computations as well as for 



those of the energies and electromagnetic forces, whose details were given in our previous 



works 



4^. 



Then we shall analyze in Sections 3 and 4 the effects of exciting a WGM in a Si cylinder 
and a LSP in an Ag one, respectively, on the optical forces exerted on them. This is carried 
out on comparing the effect of PNJ illumination with and without a subwavelength slit. In 
this way we shall show that whereas a metallic particle as small as that addressed here, 
resonantly illuminated by a nanojet, suffers an attractive vertical force, (this is in contrast 
with previous results [48(] concerning larger metallic particles), which makes the nanojet 
forming cylinder a photonic tractor, the additional presence of a subwavelength aperture 
may change the sign of this force thus transforming it to repulsive, depending on the lateral 
position of the cylinder. 

In addition, if we previously showed jl^ that the presence of the slit enhances the gradient 
forces by two orders of magnitude compared with the direct illumination of the probe particle 
by a Gaussian beam of a conventional optical tweezer, we now observe that illuminating the 
slit with a nanojet increases the magnitude of the optical forces on the particle by a factor 
between 3 and 10 compared with illumination of the slit with a Gaussian beam. 



NUMERICAL CALCULATIONS 

The calculation of the electromagnetic fields interacting with particles placed at each side 
of an aperture in a metal slab is now discussed. 2D numerical simulations are done by means 
of a finite element method (FE) (FEMLAB 3.0a of COMSOL, http://www.comsol.com). 
Aside polarization effects, the main features of the physical process of coupling, resonance 



excitation and nanojet focalization are analogous to those in 3D; this was already ana" 



for nanojets in 



14. 



38|, 



yzed 



56|. 



40l-l43| . Technical details on the calculations have been given in 
Although other materials can be chosen. Si and Ag have been used to simulate the 
probe particles, performing the analysis on the grounds of their rich Mie resonance spectra 
in the studied regimes, (near IR and near UV, respectively). The dielectric material for 
the nano jet-focusing particle is Si02 — glass |57|, l58|. The reason of working with either 
dielectric or metallic particles is to see the effects of their Mie resonances on the induced 
electromagnetic forces, taking into account the influence of the presence of a subwavelength 
slit illuminated by the nanojet. As in previous works, and because of the same reasons 



14l . l44l . l59| . the latter is considered in a thick Al slab. Noble metals like Au or may be 
used for thinner slabs. 

Following the scheme shown in Fig. 1, an incident beam, linearly p-polarized, namely 
with its magnetic vector perpendicular to the XY-plane is launched upwards, propagating 
in the Y direction. This choice is because p-polarization, in contrast with s-polarization, 
produces homogeneous eigenmodes, i. e. those which lead to extraordinary transmission 



of the subwavelength slit [381 . l44l . I60l-l62[| . The light beam incides on a dielectric cylinder 
of radius r^j, which concentrates a nanojet near its back surface. The nanojet illuminates 
either the probe cylinder of radius r when the slab in Fig. 1 does not exist; or, when the 
latter is present, this illumination is done on transmission of the nanojet through the slit, as 
seen in Fig. 1. The nomenclature followed to classify both localized surface plasmons (LSP) 
and whispering gallery modes (WGM) of the probe cylinders, as well as the resonances of 
the slit is: subscripts (i, j)^ i and j standing for their angular i-th and radial j-th orders, 
respectively. In the case of the supertransmission resonances of the slit alone the subscripts 
(u^ v) will be used, u and v standing for their longitudinal u-th and transversal v-th orders. 

In all cases the beam profile at frequency u is Gaussian, with magnetic field distribution 
at its focal plane: Hz(x,y) = \'}izo\exp{—x'^ /2a'^)exp{i{{27T/ X)y — ujt))^ |H^o| denoting the 
modulus amplitude, being the half width at half maximum (HWHM) of the beam, 

and A representing its wavelength. In this way the beam has intensity {l/7r)mW/ fim^ and 
HWHM = 2r^j. The geometrical parameters of the cylinders and of the slit have been 
adjusted such that at the given illuminating wavelength A, there is focalization of nanojets 
as well as excitation of morphology-dependent resonances of both the slit and the probe 
particle. 

The physical quantities shown in the images of the configurations under study are the 
magnetic field H^(x,^), the electric field E(x,^) and the time-averaged energy fiow < 
S(x,^) >. The curves presenting the nanojet intensity focused near the back surface of 
the Si02 — glass cylinder, are obtained at the point (x^y) where | < S(x,^) > | is maxi- 
mum. 

The light concentration either in or on the probe cylinder according to whether this is 
dielectric or metallic, is evaluated by integrating | < S(x,^) > | in a circle which coincides 
with the circle cross section of the probe cylinder of radius r when this latter is dielectric, 
or in an annulus surrounding it of radii r and rj^ when the probe cylinder is metallic, see 



FIG. 1. Schematic illustration of the geometry for the calculation of the slit transmittance and 
optical force calculations: An incident p-polarized plane wave with a Gaussian profile, (amplitude 
A = {l/7r)mW/ fim^, width of its intensity a = ^nj)? magnetic vector and time- averaged 
Poynting vector < > impinges a Si02 — glass cylinder of radius r^j which focuses a nanojet 
near its back surface, whose intensity is calculated at the point where | < > | is maximum. 
This nanojet illuminates, either directly or by interposing a subwavelength aperture of widht d, a 
cylinder of radius r. The aperture is practiced in an slab of width D and thickness h. The width 
D is also that of the horizontal side of the simulation window, whose boundary conditions are set 
to "low reflection", except in those segments which coincide with the exterior limits of the slab, 
where "perfect conductor" conditions are set. The intensity focused by the nanojet on the cylinder 
of radius r, either directly or on transmission through the slit, is evaluated as follows: For dielectric 
(Si) particles, the time averaged energy flow norm | < > | is integrated inside the cylinder circle 
cross section of radius r. For metallic {Ag) cylinders, one determines | < > | integrated in an 
annulus of exterior and interior radii: rg and r, respectively. As a beam of unit amplitude is always 
launched on the system, all transmission graphics presented in this work are normalized to this 
value. The circumference S of radius rg, is also used to calculate the electromagnetic forces on the 
probe particle (see also [14]), no matter whether the cylinder is dielectric or metallic. When there 
is no upper probe cylinder, the transmittance is obtained on integration of | < Sy{x^y) > | in a 
circle in vacuum which coincides with that of the probe cylinder cross section. 



Fig. 1. This stems from the fact that if the particle is dielectric the intensity transmitted 
by the slit, that couples to the particle WGM, is concentrated inside the cylinder; whereas 
when this nanoparticle is metallic, this transmitted intensity, coupled to a LSP, remains 
on the particle surface. On the other hand, the transmittance of the slab when the probe 
particle is not present is calculated on integration of | < S{x^y) > | in a circle of radius 
r in void space which coincides with that of the probe cylinder cross section. In all cases 
these intensities are normalized to the maximum intensity of the incident Gaussian beam 
I < Smax > I = {l/7r)mW/fim^. 



The time-averaged force on the probe cylinder is calculated by employing the Maxwell 
stress tensor (MST) [gO]: 



where the surface of integration S surrounds the particle as seen in Fig. 1 and n represents 
the outward unit normal. In our 2D geometry, S is the circumference of radius Tg, (see Fig. 



In Eq. (1), E,H and E*,H* are the values of the fields and their complex conjugates, 
6 and /i denotes the electric permittivity and magnetic permeability of the surrounding 
medium, which in this work will be assumed to be vacuum. 

The calculation with FEMLAB of the complex values E(r) and H(r) and of the real 
physical fields: E^(r,t) = Re[E{r) exp{—iujt)] and H^(r,t) = i?e[H(r) exp(— icjt)], are not 
straightforward. The details of the procedure have been given in [141] 

EXTRAORDINARY TRANSMISSION OF A SLIT IN PRESENCE OF A NANO- 
JET. EXCITATION OF A WHISPERING GALLERY MODE 

Excitation of a WGM in a cylinder by nanojet focalization. Effects of coupling by 
supertransmission 

In order to interpret the optical forces on the probe particles, we first discuss the fields 
resulting from the interaction. We begin dealing with the situation in which the slit is 




[e/2Re{{E • n)E*} - e/4(E • E*)n + /x/2i?e{(H • n)H*} 



E 



///4(H . H*)n]rf5, 



(1) 



1). e = /x= 1. 
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illuminated by a nanojet in the infrared region. The probe cylinder on the slit exit being 
dielectric. This makes it possible two processes of supertransmission enhancement: one due 
to the nanojet focused near the back surface of the dielectric cylinder of radius r^j at the 
entrance of the aperture, and that resulting from the coupling between the light exiting the 
slit and the WGM excited in the probe cylinder of radius r, (see Fig. 1). 




Wavelength (nm) 



(C) Max: 2.233W/m2(x10^) 




FIG. 2. (a) Time-averaged energy flow norm | < S(r) > |, transmitted by the slit against illumina- 
tion wavelength of: a slit alone (width d = 428. 8nm) practiced in an Alslsh (width D = 19.920/xm, 
thickness h = 857. 6nm) (black squared curve, values in right vertical axis); and of the same slit illu- 
minated by a nanojet focused near the back surface of a Si02 — glass cylinder of radius r^j = Sfim 
placed at the entrance of this slit (green triangle curve, values in left vertical axis). The red circle 
curve (values in left vertical axis) stands for | < S(r) > | concentrated by the nanojet at its max- 
imum intensity point when the slab is absent. The calculations are done as explained in Fig. 1. 
(b) Detail of the electric field |E(r)| magnitude transmission (colors) and vector (arrows) of the 
nanojet focused near the back surface of the Si02 — glass cylinder at A = 1175nm, without slit 
nor Si probe particle. This corresponds to a peak of the red circle line of (a), (c) The same detail 
as in (b) but of the time- averaged energy flow < S(r) >, (norm in colors and vector in arrows) . 

We first make a study of the fields and energy without the probe Si. Fig. 2(a) shows 
the black squared and green triangled curves which stand for the average energy flow, either 
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when the Gaussian beam directly illuminates the Al slab without the presence of the nanojet, 
(peak at A = 1175nm), or when the incident light is focused as a nanojet near the back 
surface of a Si02 — glass cylinder of radius Tnj — 3/im placed at the slit entrance. The 
intensity at that point of the nanojet where it is maximum, without slab, is plotted as the 
red circle curve. The presence of the nanojet enhances the transmittance peak of the slit by 
a factor above 4 (at A = 1155nm), but erases the supertransmitted energy peak profile of 
the slit alone (compare the green triangle and the black square lines), manifesting the non 
resonant nature of the PNJ focusing process (see the red circle hne). A detail of the nanojet 
alone, i.e. without slit nor Si cylinder, near the back surface of the cylinder is displayed 
both by its electric field E(r) distribution [see Fig. 2(b)] and by that of its energy fiow 
< S(r) > [Fig. 2(c)], both maps clearly show the focalization process and the interesting 
vortex of electric field lines about the nanojet focus. 

Next, the Si probe cylinder is placed above the Si02 — glass without the presence of 
the slit. The light in the Si cylinder now concentrates in the form of a WGH2^i resonance, 
as shown by the dipolar map of the magnetic field norm |H(r)| in Fig. 3(a), r = (x,y). 
Consequently, the electric field E(r) distribution acquires a vortex pattern near the probe 
cylinder [cf. Fig. 3(b)] and, under this incident p-polarization, the light intensity from the 
nanojet is strongly localized within this Si cylinder [see Fig. 3(c)]. 

Figure 4(a) shows the average energy flow for all the conflgurations studied in this paper. 
The intensity of the nanojet localized in the Si cylinder of radius Tnj — 3fim^ without slit, is 
plotted in the blue curve with triangles down. The peaks of this curve coincide with those of 
the intensity peak point of the nanojet alone, [see the red curve with circles in Fig. 2(a)], but 
now enhanced by the resonant concentration due to the WGH2^i resonance at A = 1175nm. 
As seen, this nanojet also enhances the energy asociated to the cylinder resonance by a factor 
of 10, [cf. the blue curve with triangles down in Fig. 4(a)], on comparing with the case of 
the resonant Si cylinder directly illuminated by the Gaussian beam without nanojet nor slit, 
(whose intensity integrated in the Si cylinder cross-section has a maximum at A = 1195nm; 
not shown here for brevity). 

Concerning the nanojet focused in the Al slit with the Si cylinder at its exit, the trans- 
mitted intensity inside that cylinder against wavelength takes the form of the pink curve 
with rotated triangles of Fig. 4(a). This follows the trend of the intensity in that cyhnder 
alone, (not shown here), modulated by that of the nanojet-slit system, as a result of the 
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(a) Max: 2.71 9A/m(x10^) (b) Max: 2.411 V/m(x10^) 




FIG. 3. (a) Detail of the map of the magnetic field norm |H(r)|, r = (x, y), in a resonant Si cylinder 
of radius r — 200nm illuminated by the nanojet focused at 350nm near the back surface of the 
Si02 — glass cylinder of radius r^j = 3/xm. The Si cylinder is placed at the maximum intensity 
point of the nanojet. There is energy concentration as a WGH2^i mode of the Si cylinder, (b) 
Detail of the electric field E(r) magnitude and vector, (colors and arrows, respectively), (c) Detail 
of time-averaged energy flow < S(r) >, (norm in colors and vector in arrows). The conditions of 
illumination in these images are those of Figs. 2(b) - (c) and at the same wavelength. 

effect of nanojet focalization plus WGM resonance excitation. This shows a modest increase 
of transmitted intensity inside the Si particle (cf. the pink curve with rotated triangles) 
comparing with that of the nanojet large cylinder plus slit, without the Si particle, (cf. the 
green curve with triangles up). Nevertheless, the enhancement factor is now above 4 (at 
A = 1160nm) with respect to the transmitted intensity from the slit alone, [cf. black curve 
with squares in Fig. 2(a)]. 

A detail of the spatial distributions of E and < S > at the optimum wavelength for 
largest transmitted intensity of the system constituted by the large Si02 — glass cylinder, 
the subwavelength slit and the resonant Si probe cylinder, can be seen in Figs. 4(b) and 
4(c), respectively, showing the vortex of the electric vector in the Si cylinder around the 
point of maximum intensity concentration. Also, the large charge localization in the corners 
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(b) Max: 5.00V/m(x10^) 




FIG. 4. (a) Transmitted time-averaged energy flow norm | < S(r) > | versus illumination wave- 
length of the following configurations: A nanojet focused by the Si02 — glass cylinder of radius 
Tnj — tangent to the entrance plane of the slit practiced in the Al slab studied in Fig. 2(a), 

(green curve with up-triangles) . A probe Si cylinder is added, tangent to the exit plane of the slit 
illuminated by the nanojet (pink curve with rotated triangles). The blue curve with down-triangles 
stands for | < S(r) > | concentrated in the resonant Si cylinder, placed above the nano jet-focusing 
Si02 — glass one, when there is no slab between them. All these curves are obtained as explained 
in Fig. 1. (b) Detail of the electric field E(r) magnitude and vector (colors and arrows) in the Si 
cylinder at the exit plane of the slit illuminated by the nanojet, (see in (a) the pink curve with 
rotated triangles). The illumination of the nanojet forming large Si02 — glass cylinder below the 
slit is with a Gaussian beam at A = llGOnm. This cylinder is not shown in the figure, (c) The 
same as in (b) for < S(r) >, norm in colors and vector in arrows. 



of the slit is clearly observed. This will be seen to have important consequences for the 
optical forces on the Si cylinder. 
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Electromagnetic forces on a dielectric cylinder with a whispering gallery mode 
excited by a nanojet, either directly or on transmission through a subwavelength slit. 



Next, we study the force that the electromagnetic fields of the nanojet transmitted 
through the slit exert on the probe Si cylinder placed at the slit exit. 
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FIG. 5. (a) X- and Y-components of the total electromagnetic force on the Si cylinder as it 
laterally moves to the right from the point of maximum intensity of the nanojet focused by the 
Si02 — glass cylinder of radius rnj = 3fim. This point is at 350nm from the back surface of this 
Si02 — glass cylinder. No slit is present, (b) The same quantities when the Si cylinder separates 
vertically from the same point as in (a). In both graphics, black squares and green up-triangles 
curves stand for when the Si cylinder is non-resonant (A = 1280nm) and when it is resonant 
(A = 1175nm), respectively. Red circle and blue down-triangles represent, in the same order, F^. 



Figure 5(a) shows the variation in the X- and Y-components of this force when the 
particle moves horizontally from the peak point of the nanojet intensity map. In this case 
we assume no slit present. The particle, illuminated out of resonance, is gradually and 
increasingly attracted along the lateral direction by the gradient of intensity of the nanojet 
spatial distribution, and less and less pushed vertically, (cf. black squared and red circled 
curves, respectively). On the other hand, when the WGH2^i is excited, the particle is now 
slightly more and more laterally repelled from the peak intensity point of the nanojet map, 
(see the green up-triangle curve) and weakly, although increasingly and progressively, pushed 
in the vertical direction (cf. blue-down triangle curve). 

Figure 5(b) presents the same quantities as those plotted in Fig. 5(a), but with the Si 
cylinder vertically moving upwards from the maximum intensity point of the nanojet spatial 
distribution. The absence or presence of the resonance in this particle is shown by the black 
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square and green up-triangle curves, respectively. Obviously, due to the symmetry of the 
field intensity map around this axis, no horizontal force is exerted. In the vertical direction 
the particle is nevertheless fairly strongly pushed upwards, the force magnitude following an 
oscillatory behavior due to the interference pattern between the field exiting the slit and that 
refiected down by the Si cylinder. This push is weak if the WGM resonance is excited. In 
this latter case, under this p-polarization illumination most of the near field is concentrated 
inside the particle, with little intensity outside. This leads to such weak optical forces on 
the resonant particle with the p-polarized WGM in both Figs. 5(a) and 5(b). 

It is worth remarking that these forces and induced by the nanojet on the Si particle 
are about 1000 and 10 times larger, respectively, than those from direct illumination of this 
probe cylinder by the Gaussian beam like in a conventional optical tweezer, (compare the 
above Figs. 5(a) and 5(b) with Figs. 5(a) and 5(b) of fl^). 
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FIG. 6. (a) X-component F^ of the total electromagnetic force exerted on the Si cylinder placed 
on the exit of the slit practiced in the Al slab illuminated by a nanojet that arises near the back 
surface of the Si02 — glass cylinder of radius r^j = 3fim below the slit, See Fig.l. The Si cylinder 
is tangent to the exit plane of the slit, and laterally moves to the right from the aperture axis, 
(b) Y-component F^ of the optical force when the Si cylinder moves as in (a). In both graphics, 
black square and red circle curves stand for the cases in which the probe cylinder is non-resonant 
(A = 1280nm) and resonant (A = llGOnm), respectively. 



On the other hand. Figs. 6(a) and 6(b) show the X- and Y-components of the optical 
force exerted on the Si cylinder, placed on the exit of the subwavelength slit, by the electro- 
magnetic field exiting it. A Si02 — glass cylinder of radius r^j = 3/xm is below the slit and 
forms a nanojet near its back surface. The particle laterally moves to the right approaching 
the slit edge. Each force Cartesian component is plotted for the Si cylinder being either out 
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or in resonance. Now this particle is more and more horizontally attracted by the corners of 
the slit exit as it approaches them, independently of the absence or presence of the resonance 
[see the black square and red circle curves, respectively, in Fig. 6(a)], the horizontal force 
being stronger in resonance. Notwithstanding, the slit exerts a pushing vertical force on the 
cylinder which decreases as it moves to the slit corners, [cf. black squared and red circled 
curves, respectively, in Fig. 6(b)]. Then, the presence of the resonant WGM2^i in the Si 
particle causes the vertical push from the slit to decrease and the horizontal attraction to 
the slit corners to increase. The cause of this is the same as discussed in connection with 
Figs. 5(a) and 5(b) for p-polarized illumination. However, as regards the horizontal force 
near the slit edge, the gradient force created around this edge by its concentration of charge, 
[see Figs. 4(b) and 4(c)], is responsible for the larger value of in resonance near this 
edge, as seen in Fig. 6(a). Thus, solely the field intensity gradient near the exit corners of 
the slit, (remember that in resonance under p-polarization the Si cylinder concentrates the 
intensity inside, leaving very little scattered energy in its surroundings), is the major factor 
concerning the interaction between the aperture and the resonant particle. 

These forces induced by the nanojet on the Si particle through the slit are about 3 and 
6 times larger than those from direct illumination of the slit by the Gaussian beam, i.e. 
without nanojet, (compare the above Figs. 6(a) and 6(b) with Figs. 5(a) and 5(b) of [14]). 



EXTRAORDINARY TRANSMISSION OF A SUBWAVELENGTH SLIT ILLUMI- 
NATED BY A NANOJET. EXCITATION OF A LOCALIZED SURFACE PLAS- 
MON. 

Excitation of a localized surface plasmon in a cylinder by nanojet focalization. 
Effects of coupling by supertransmission 

We now deal with a subwavelength slit illuminated by a nanojet in the ultraviolet region 
when a metallic cylinder is at its exit. This produces two interacting processes of extraor- 
dinary transmission, or supertransmission, enhancement: one due to the coupling of the 
nanojet focused in the aperture with the excitation of its propagating eigenmodes, and that 
arising from the coupling of the field transmitted through the slit with the localized surface 

nn 

plasmon (LSP) excited on the metallic cylinder, [IJ, In order to analyze the optical 
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force on the metallic particle, we first briefly discuss the distribution of flelds and energy in 
this system. 
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FIG. 7. (a) Time-averaged energy flow norm | < S(r) > | versus illumination wavelength, trans- 
mitted at the exit of a slit (width d = 109. 7nm) practiced in an Al slab (width D = 5.096/im, 
thickness h = 219. 4nm) illuminated by a nanojet. This nanojet is focused by a Si02 — glass 
cylinder of radius Vnj = 1.3/xm, lOnm separated from the entrance plane of the slit, (b) No slit is 
present: | < S(r) > | against illumination wavelength, concentrated on the surface of an Ag cylinder 
of radius r — 30nm placed at 200nm from the back surface of the focusing Si02 — glass cylinder, 
(c) The same quantity when the slit is present and the Ag cylinder, tangent to the exit plane of the 
aperture, scatters the light exiting the slit which is illuminated by the nanojet. These calculations 
are done as explained in Fig. 1. 



Now a Si02 — glass cylinder of radius Tnj — 1.3/im is placed below an Al slab and focuses, 
near its back surface, a nanojet which illuminates the slit. The calculations are done as 
described in Fig. 1. The rather smooth transmitted intensity with peak at A = 339. 7nm 
that would be obtained if this slit were illuminated by the Gaussian beam without any 
nanojet, normalized to the maximum intensity (which was 0.75mW/ fim^)^ becomes, when 
the nanojet illuminates the slit, an oscillatory line, typical of nanojets, as seen in Fig. 7(a). 
Nevertheless, this conflguration may render a slit-nanojet transmittance enhancement factor 
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of about 4. Figure 7(b) represents the effect of light concentration in an Ag cyhnder of radius 
r = 30nm placed at 200nm from the back surface of the focusing Si02 — glass cylinder when 
the slit is suppressed. The excited TM2^i eigenmode energy concentrated on the surface of 
the Ag cylinder is enhanced at A = 339. 7nm by a factor about 33 due to the presence of 
the nanojet, [see Fig. 7(b)]. This effect is somewhat similar to that on the above studied Si 
particle, although now with a much stronger enhancement factor. 

On the other hand, when the slit is also present as shown in Fig. 7(c), enhancement 
factors between 8 and 47 are obtained on ecitation of the TM2,i LSP in the Ag cylinder. 
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FIG. 8. (a) Detail of the electric field spatial distribution E(r), r = {x,y), (magnitude and 
vector: colors and arrows, respectively), near the Ag cylinder of radius r = 30nm placed at 200nm 
from the back surface of the focusing Si02 — glass cylinder of radius r^j = l.Sfim. No slit is 
present, (b) The same detail but now in time- averaged energy flow < S(r) >, (norm in colors 
and vector in arrows), (c) Detail of E(r) (magnitude and vector: colors and arrows, respectively), 
when the slit is also present. The configuration corresponds to the largest peak in of Fig. 7(c). 
|E(r)|^^^ ~ 2.513 • lO^V/m. (d) The same detail as in (c), now in < S(r) > (norm in colors and 
vector in arrows). | < S(r) > |^^^ = 9.265 • lO^^W/m^ . In all images the illumination wavelength 
is A = 339. 7nm. 



As illustrations. Figs. 8(a) and 8(b) present details of the effect of the nanojet near the 
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Ag cylinder. As seen, there is a skin depth effect at the illuminating wavelength due to the 
small size of this Ag cylinder. This case corresponds to the largest peak in Fig. 7(b). The 
electric field E(r) and the time- averaged energy fiow < S(r) >, respectively, are also shown. 
The nanojet is almost completely concentrated in the small metallic cylinder. Figs. 8(c) 
and 8(d) are details of E(r) and < S(r) >, respectively, for the configuration studied in Fig. 
7(c) i. e. that of the nanojet-LSP coupling through the slit. The wavelength A = 339. 7nm 
corresponds to the largest peak in Fig. 7(c), which practically remains the same as that of 
Fig. 7(b) when there is no slit. The hotspots in the slit edge, the nanojet focalization as 
well as the LSP excitation mechanisms are shown. 



Electromagnetic forces on a metallic cylinder with a localized surface plasmon ex- 
cited by a nanojet, either directly or on transmission through a subwavelength slit. 

We now are in position to analyze the optical forces induced by the nanojet and the slit 
on the metallic Ag particle with excitation LSPs. 



(a) 



(b) 
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FIG. 9. (a) X- and Y-components of the total electromagnetic force F^, calculated as explained 
in Fig. 1 and Eq. (1), due to the nanojet focused by the Si02 — glass cylinder, acting on the 
Ag cylinder of radius r = 30nm as this laterally moves to the right from the point of maximum 
intensity of the nanojet. (This point is on the vertical symmetry axis at 200nm from the back 
surface of the Si02 — glass cylinder). No slit is present, (b) The same when the Ag cylinder 
vertically moves upwards from the point of maximum intensity of the nanojet. In both figures, 
black squared and triangle up green curves stand for F^, both when the Ag cylinder is non-resonant 
(A = 500. Onm) and when it is resonant (A = 339. 7nm), respectively. Red circle and triangle down 
blue lines represent F^ in the same order and for the same cases as before. 
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To this end, the excitation of the Ag cyhnder TM2^i LSP eigenmode by focusing a nanojet 
on it, is studied. Figure 9(a) shows the behavior of the X- and Y-components of the force on 
this metaUic particle as it progressively moves to the right in absence of the slit. This force 
is due to only the field of the nanojet scattered by the Ag particle. When the particle is not 
resonant, the forces along the X and Y directions are very weak, of the order of 10~'^fN/nm 
and 10~^fN/nm for the horizontal force (black curve with squares) and for the vertical force 
(red curve with circles), respectively. In resonance the horizontal force is still very weak, 
but the vertical component is strong and atractive. On the other hand, as the Ag cylinder 
moves upwards along the nanojet axis , as shown in Fig. 9(b), no transversal force appears 
(black squared and green up-triangled curves, respectively) as it should, but differences in 
the vertical force arise depending on the absence or the presence of a LSP resonance. This 
vertical force is very small in absence of LSP excitation: about 10~^fN/nm^ whereas it 
is attractive and oscillating with distance, and of the order of 10^ fN/nm when the LSP 
resonance of the Ag cylinder is excited, this oscillation is due to the interference pattern 
between the nanojet inciding on the metallic particle and the field that the latter scatters 
backward. 

Notice the remarkable pulling nature of this vertical force, which converts the Si02— glass 
cylinder in a photonic tractor 63- 65| for certain plasmonic small particles. This contrasts 
with the repulsive obtained in [48] on an Ag particle in resonance. However, if one 
doubles its radius making it of 60nm, then this force becomes repulsive due to its larger 
scattering cross section. We do not show this here for brevity. 

Finally, the presence of the slit illuminated by the nanojet is analized concerning the 
optical forces acting on the Ag cylinder. Once again, the excitation of a LSP resonance 
largely aflFects these forces, but due to p-polarized illumination, it has consequences opposite 
to those of the WGM excitation in a dielectric particle, (cf. Figs. 10(a) and 10(b) for 
the X- and Y-components of this optical force, respectively). Black square and red circle 
curves in Fig. 10(a) stand for the cases of absence and presence of LSP resonance in the 
Ag cylinder, respectively. They show that the plasmonic resonance enhances the lateral 
force, which progressively attracts the particle towards the corners as it approaches it in its 
displacement. The vertical force, plotted with the same color code in Fig. 10(b), behaves 
in a more complex way. This being slightly attractive when the Ag cylinder is non resonant 
(values of order 10~^ fN/nm). However, in presence of resonance this vertical force is, 
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FIG. 10. (a) X-component of the total electromagnetic force exerted on the Ag cylinder of radius 
r = 30nm, calculated as explained in Fig. 1 and Eq. (1). The slit is practiced in the Al slab, 
being illuminated by the nanojet which arises from the back surface of the Si02 — glass cylinder 
placed below. The Ag cylinder is tangent to the exit plane of the slit, and laterally moves to the 
right from the aperture axis, (b) Y-component of the optical force when the Ag cylinder moves 
as in (a). In both graphics, black square and red circle curves stand for the cases in which the 
Ag cylinder, tangent to the exit line of the aperture, is non-resonant (A = 1280nm) and resonant 
(A = llGOnm), respectively. The green up triangles and blue down triangles curves represent, in 
the same order as before, the same cases when the Ag cylinder is at the exit slit, 20nm inwards. 

though, repulsive in the central region of the slit exit, becoming attractive near its corners 
(red circled curve). These tendencies can be increased by locating the Ag cylinder inwards 
the slit exit, as noticed by the green up-triangles and blue down-triangles curves of both 
Figs. 10(a) and 10(b). [Values of order 10~^fN/nm for the green up triangles curve (non 
resonant cylinder) in Fig. 10(b)]. 

In any case, the forces induced by the nanojet on the Ag particle through the slit are 
about 10 times larger than those from direct illumination of the slit by the Gaussian beam 
without nanojet, (compare the above Figs. 10(a) and 10(b) with Figs. 11(a) and 11(b) of 

DISCUSSION AND CONCLUSIONS 

We have made a study on the photonic forces on nanoparticles near a subwavelength 
slit iUuminated by a nanojet. We have performed comparisons of the force on the particles 
either in or out of their Mie resonance excitation. The study has been done in 2D, so that 
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these objects are cylinders and the transmission modes of the sht are excited by p-polarized 
iUumination. However, we beheve that these results also hold for spheres in 3D, aside of 
polarization effects. 

We have proven that whereas the morphological resonance of the slit, causing supertrans- 
mission, enhances the fields surrounding the probe cylinders and hence the optical forces 
exerted upon them, the excitation of the probe particle resonance has a different effect on 
these forces, according to whether the cylinder is dielectric or metallic. 

The photonic force on the probe particle dependens on the field surrounding it. Thus 
WGMs in dielectric nanocylinders under p-polarization, (which is the one at which super- 
transmission exists in 2D), have no much effects on the optical force. Its strength may 
weaken due to the low intensity distribution in the near field region of this particle. In 
this connection, it must be remarked that the electromagnetic force on dielectric cylinders 



50|; then the WGM extends to the 



is however enhanced under s-polarization illumination 
near field outside the cylinder. 

On the other hand, LSPs in metallic nanocylinders give rise to localized field energy 
enhancement on the cylinder surface; this strengthes the optical force. In contrast with 
illumination by only a PNJ, the presence of a subwavelength slit enriches the effect of the 
electromagnetic forces. Whereas such a a small metallic particle resonantly illuminated by 
a nanojet, suffers an attractive vertical force which makes the Si02 — glass nanojet forming 
cylinder to behave as a photonic tractor, the additional presence of a subwavelength aperture 
may change the sign of this force thus transforming it to repulsive, depending on the lateral 
position of the cylinder. 

For any cylinder, whether dielectric or metallic, we observe that illuminating the slit with 
a nanojet increases the magnitude of the optical forces on the particle by a factor between 
3 and 10 compared with illumination of the slit with a Gausian beam. On the other hand, 
the gradient horizontal force induced by a nanojet in a dielectric particle is a few orders of 
magnitude larger than that of a Gaussian beam of a conventional optical tweezer. 

These results should stimulate experiments with PNJs on slits since they are the result 
of focusing by microparticles, and their highly localized intensity at the subwavelength scale 
present an ideal means to spatially control the mechanical action on nanoobjects. 
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